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I. INTRODUCTION 

The quest for a quantum theory of gravity involves searching for a microscopic quantum theory whose low energy limit is the 
known physics of general relativity, dynamical space-time metrics whose evolution is governed by the Einstein equations. Many 
approaches, like loop quantum gravity JT|, Causal Dynamical Triangulations (CDT) [2|, spin foams [3] and group field theory 
0, expect that the quantum theory of gravity becomes manifest at very high energy. That is, quantum analogs of geometry and 
gravitational properties such as the quantum Hilbert-Einstein action or Lorentz invariance are built into the high energy theory. 
In different theories, various such properties are present at the microscopic level: In loop quantum gravity, for example, the 
state space is given in terms of spin network basis states, understood as a discretization of space (see 0|6l[8)) and, in addition, 
embedded in 3-dimensional space. The dynamics of the micro-states is governed by the Hamiltonian constraint, obtained by 
the canonical quantization of the Hilbert-Einstein action. In spin foams or group field theory, the microscopic states similarly 
are based on simplicial graphs or complexes carrying algebraic geometric data, including the Lorentz group. Depending on the 
model, a spin foam may be based on states that are embedded complexes or combinatorial (abstract, non-embedded graphs). In 
many models, the states carry representations of the Lorentz group, while the dynamics uses the Regge action, a discrete form 
of the Hilbert-Einstein action. In CDT, the states are simplicial combinatorial objects, again encoding microscopic Lorentzian 
geometry, while the dynamics contains the Regge action. Of course, the important question is how many of these features survive 
at the low energy or continuum limit. Is the input of geometric and gravitational properties at the microscopic level necessary 
for the continuum limit to be gravity? The answer is at present unclear. CDT results indicate that while certain properties 
such as causality are important, others, including the gravitational action, may not be (9). In spin foams or group field theory, 
it is not clear that the Lorentz group representations on the microstates ensure the reappearance of the Lorentz group in the 
continuum limit. As an alternative, one could regard quantum gravity as an emergent phenomenon from the low energy theory 
of a condensed matter system. In this approach, the fundamental theory would have no microscopic degrees of freedom, no 
gravity, no elements of the Lorentz group, etc. 

In recent years, there has been a growing attention to the notion of gravity as an emergent phenomenon 1101 . From Aristotle 
to Philip Anderson, a long-standing tradition in physics asserts that "The whole is more than the sum of its parts" and that "More 
is different". The emergent approach is concerned with the study of the macroscopic properties of systems with many bodies. 
Sometimes, these properties can be tracked down to the properties of the elementary constituents. In recent years, though, there 
has been a flourishing of novel quantum systems, which show behaviors of the whole system that have no explanation in terms 
of the constituting particles, but instead of their collective behavior and interaction. When the interaction between the particles 
cannot be ignored, like in systems of strongly interacting electrons, we see many novel and beautiful properties: gauge fields 
can emerge as a collective phenomenon, strange quantum phase transitions happen, unusual forms of superconductivity and 
magnetism appear, novel orders of the matter based on topological properties of the system and featuring exotic statistics are 
found. 

In this sense, one can view the problem of quantum gravity as a problem in statistical physics or condensed matter theory: we 
know the low energy physics and are looking for the correct universality class of the microscopic quantum theory. By analogy 
to the Ising model for ferromagnetism, one can ask: What is the "Ising model" for gravity? A number of such approaches to 
the problem have recently been proposed, ranging from looking for gravity analogues in condensed matter systems fl2l[T3l . to 
condensed matter systems with emergent graviton excitations [ 14|, and spin system models for emergent geometry |[T6l[T7ll . One 
can also view approaches such as CDT |2 1 or matrix models [ 18 1 in this light. In comparison to the quantum gravity approaches 
mentioned above, in the systems we will be investigating there is no straightforward geometric content to the microscopic 
graph states. We simply use a dynamical network of quantum relations to describe a world without geometry, where locality is 
determined by the presence of absence of quantum interactions. It is perhaps best to think of this system as a quantum information 
processor, with information more primary than geometry, and geometry being the set of properties such as geometric symmetries 
that we only expect to find dynamically emergent at low energy. 

An important issue in this direction of research in quantum gravity is the dynamical nature of geometry in general relativity. 
Normally, methods in condensed matter theory use a fixed background, for example, a spin system on a fixed lattice. The lattice 
determines the locality of the interactions and hence is a discretization of geometry. One can then worry that this direction of 
research is limited to fixed background geometries. This is addressed in some, but not all, of the above-mentioned approaches. 
For example, CDT is based on path integral dynamics on an ensemble of all lattices (each lattice is a regularization of a Lorenzian 
geometry), thus providing a proper non-perturbative approach to the problem. Elsewhere, ideas from quantum information theory 
have been introduced to deal with this problem ifTTl |2T1 l22l . 

An alternative direction, and the one we are pursuing in the present article, is to make the lattice itself dynamical. Two of 
the present authors proposed such a model for the emergence of geometry in [16|. The basic idea was to promote the lattice 
links to dynamical quantum degrees of freedom and construct a Hamiltonian such that, at low energy, the system "freezes" in 
a configuration with recognizable geometric symmetries, interpreted as the geometric phase of the model. The present work 
revisits the same idea but in a different model with two central properties: 

• The model is a spin system on a dynamical lattice. 



3 



• There are lattice and matter degrees of freedom. The lattice interacts with the matter: matter tells geometry how to curve 
and geometry tells matter where to go. 

The starting point for the implementation of the above is considerations of locality. Normally, locality is specified by the metric 
g^y on a manifold A4. Dynamics of matter on (A4,g^) is given by a Lagrangian which we call local if the interaction terms 
are between systems local according to g^ u . A Lagrangian with non-local interaction terms is typically considered unphysical. 
That is, the matter dynamics is made to match the given space-time geometry. We will do the reverse and define geometry via 
the dynamics of the matter. Our principle is that if particles i and j interact, they must be adjacent. This is a dynamical notion 
of adjacency in two ways: it is inferred from the dynamics and, being a quantum degree of freedom, it changes dynamically in 
time. This amounts to a spin system on a dynamical lattice and to interaction of matter with geometry. 

To summarize, we present a toy model for the emergence of locality from the dynamics of a quantum many-body system. 
No notion of space is presupposed. Extension, separateness, distance, and all the spatial notions are emergent from the more 
fundamental notion of interaction. The locality of interactions is now a consequence of this approach and not a principle. We 
will promote the interaction terms between two systems to quantum degrees of freedom, so that the structure of interactions 
itself becomes a dynamical variable. This makes possible the interaction and even entanglement between matter and geometry. 

This toy model is also a condensed matter system in which the pattern of interaction itself is a quantum degree of freedom 
instead of being a fixed graph. It can be regarded as a Hubbard model where the strength of the hopping emerges as the mean 
field value for other quantum degrees of freedom. We show a numerical simulation of the quantum system and results on the 
asymptotic behavior of the classical system. The numerical simulation is mainly concerned with the entanglement dynamics of 
the system and the issue of its thermalization as a closed system. A closed system can thermalize in the sense that the partial 
system shows some typicality, or some relevant observables reach a steady or almost steady value for long times. The issue of 
thermalization for closed quantum system and the foundations of quantum statistical mechanics gained recently novel interest 
with the understanding that the role of entanglement plays in it |38|. The behavior of out of equilibrium quantum system under 
sudden quench, and the approach to equilibrium has been recently the object of study to gain insight in novel and exotic quantum 
phases like topologically ordered states. 

From the point of view of Quantum Gravity, the interesting question is whether such a system can capture aspects of the 
dynamics encoded by the Einstein equations. We start investigating in this direction by studying an analogue of a trapped 
surface that may describe, in more complete models, black hole physics. We discuss physical consequences of the entanglement 
between matter and geometry. 

The model presented here is very basic and we do not expect it to yield a realistic description of gravitational phenomena. 
What we would like to show is that such a model can have an emergent, quantum-mechanical notion of geometry (even if not 
smooth), that locality is derivative from dynamics, and the extent to which such a simple model may capture aspects of the 
Einstein equations, is left to future work. 



n. THE MODEL WITH HOPPING BOSONS ON A DYNAMICAL LATTICE 

A. Promoting the edges of the lattice to a quantum degree of freedom 

We start with the primitive notion of a set of N distinguishable physical systems. We assume a quantum mechanical de- 
scription of such physical systems, given by the set {Hi, Hi} of the Hilbert spaces Hi and Hamiltonians Hi of the systems 
i = 1, N. This presumes it makes sense to talk of the time evolution of some observable with support in Hi without making 
any reference to space. 

We choose Hi to be the Hilbert space of a harmonic oscillator. We denote its creation and destruction operators by b\,bi, 
respectively, satisfying the usual bosonic relations. Our N physical systems then are N bosonic particles and the total Hilbert 
space for the bosons is given by 

N 

Hbosons — ^^Hi- (1) 
i=l 

If the harmonic oscillators are not interacting, the total Hamiltonian is trivial: 

N 

H v = Y,Hi = -Y,»Ak- (2) 

i— 1 i 

If, instead, the harmonic oscillators are interacting, we need to specify which is interacting with which. Let us call 1 the set of 
the pairs of oscillators e = (i, j) that are interacting. Then the Hamiltonian would read as 



H = £ Hi + ^ 



(3) 
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where h e is a Hermitian operator on %i ® Wj representing the interaction between the system i and the system j. 

We wish to describe space as the system of relations among the physical systems labeled by i. In a discrete setup like ours, a 
commonly used primitive notion of the spatial configuration of N systems can be provided by an adjacency matrix A, the N x N 
symmetric matrix defined as follows: 

. J 1 if i and j are adjacent 

lj ~\0 otherwise. ' (4) 

The matrix A is associated to a graph on N vertices whose edges are specified by its the nonzero entries. Now, it is clear that 
the set 1 of interacting nodes in the Hamiltonian ([3]) also defines a graph G whose vertices are the N harmonic oscillators and 
whose edges are the pairs e = (£, j) of interacting oscillators. Here I is the edge set of G. We want to promote the interactions 
- and thus the graph itself - to a quantum degree of freedom. 

To this goal, let us define Q as the set of graphs G with N vertices. They are all subgraphs of Km, the complete graph on N 
vertices, whose JV ( J ^~ 1 ) edges correspond to the (unordered) pairs e = of harmonic oscillators. To every such pair e (an 
edge of Kn) we associate a Hilbert space 7i e ~ C 2 of a spin 1/2. The total Hilbert space for the graph edges is thus 

jV(AT-l)/2 

Hgraph = ^e- (5) 

e=l 

We choose the basis in H gra ph so that to every graph g E Q corresponds a basis element in H gra ph' the basis element 
|ei . . . eM(N-i)/2) = \G) corresponds to the graph G that has all the edges e s such that e s — 1. For every edge the 
corresponding SU(2) generators will be denoted as S l = l/2a l where a 1 are the Pauli matrices. 
The total Hilbert space of the theory is 

it ^-bosons ^ i^graphi (^) 

and therefore a basis state in H has the form 

|^ = |^(6 oso „ s)) ^( 9 ra ph)) = |„ ^g, | e g^.,,) (7 ) 

2 

The first factor tells us how many bosons there are at every site i (in the Fock space representation) and the second factor tells 
us which pairs e interact. That is, the structure of interactions is now promoted to a quantum degree of freedom. A generic state 
in our theory will have the form 

i$> = E^i vI 'i b ~ ) )^i*i 9rap ' l) )' (8) 

a, 6 

with ^ a b |a a ,b| 2 = 1. In general, our quantum state describes a system in a generic superposition of energies of the harmonic 
oscillators, and of interaction terms among them. A state can thus be a quantum superposition of "interactions". For example, 
consider the systems i and j in the state 

= |10) + |io) sio),,- (9) 

V2 

This state describes the system in which there is a particle in i and no particle in j, but also there is a quantum superposition 
between i and j interacting or not. The following state, 

= |00) <E> + |11) (8> |0>»j _ (10) 
v2 

represents a different superposition, in which the particle degrees of freedom and the graph degrees of freedom are entangled. It 
is a significant feature of our model that matter can be entangled with geometry. 

An interesting interaction term is the one that describes the physical process in which a quantum in the oscillator i is destroyed 
and one in the oscillator j is created. The possibility of this dynamical process means there is an edge between i and j. Such 
dynamics is described by a Hamiltonian of the form 

H hop = -tJ2 P H ® ( b t b 3 + b 4) ^ 
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where 

is the projector on the state such that the edge is present and the spin operators are defined as = |l)(0|(.y) and 

^(if) = With this Hamiltonian, the state Eq.(|9j can be interpreted as the quantum superposition of a particle that 

may hop or not from one site to another. It is possible to design such systems in the laboratory. For instance, one can use arrays 
of Josephson junctions whose interaction is mediated by a quantum dot with two levels. 

We note that it is the dynamics of the particles described by Hf lop that gives to the degree of freedom \e) the meaning of 
geometry [? ]. The geometry at a given instance is given by the set of relations describing the dynamical potentiality for a 
hopping. Two points j, k can be "empty", that is, the oscillators j, k are in the ground state, but they can have a spatial 
relationship consisting in the fact that they can interact. For example, they can serve to have a particle to hop from i to j, then to 
k, then to I, We read out the structure of the graph from the interactions, not from the mutual positions of particles. 

In addition, H^op tells us that it takes a finite amount of time to go from i to j. If the graph is represented by a chain, it tells 
us that it takes a finite amount of time (modulo exponential decaying terms) for a particle to go from one end of the chain to 
another. This results to a "spacetime" picture (the evolution of the adjacency graph in time) with a finite lightcone structure. The 
hopping amplitude is given by t, and therefore all the bosons have the same speed. We can make the model more sophisticated 
by enlarging the Hilbert space of the links, and obtain different speeds for the bosons. Instead of considering spins 1/2, consider 
an S— level system. The local Hilbert space is therefore 

H e = span{|0),|l),...,|S-l)} (13) 

Now consider the projector onto the s— th state on the link Pffl = \a)(s\ij. We can define a new hopping term whose 

amplitude depends on the level of the local system in the following way: 

H hop = - tsP ij ] ® ( b i b 3 + bib \) (14) 

where the hopping amplitudes t s depend on the state s of the system, and to = 0. For instance, the t s can be chosen larger for 
larger s. In this way, moves through higher level links are more probable, and therefore the speed of the particles is not constant. 
In the following, we will study the model with just the two level system. 

Of course, we need a Hamiltonian also for the spatial degrees of freedom alone. The simplest choice is simply to assign some 
energy to every edge: 

Finally, we want space and matter to interact in a way that they can be converted one into another. The term 

H ex = kJ2 (S{ id) ® {b\b)) R + S+ j} ® (hb 3 ) R ) (16) 

can destroy an edge (i, j) and create R quanta at i and R quanta at j, or, vice-versa, destroy R quanta at i and R quanta at j to 
convert them into an edge. 

The terms Hunu and H ex are so simple that we will not expect them to give us any really interesting property of how regular 
geometry can emerge in such a system. This is the subject for a more refined and future work. Nevertheless, this term has 
an important meaning because the nature of the spatial degrees of freedom is completely reduced to that of the quanta of the 
oscillators: an edge is the bound state of 2R quanta. When in the edge form, the quanta cannot hop around. When unbounded, 
they can hop around under the condition that there are edges from one vertex to another. One can replace the separation of the 
fundamental degrees of freedom into bosons and graph edges with a unified set of underlying particle ones, single bosons and 
collections of 2R bound bosons. Therefore a bound state of 2R quanta in the pair (Hi,Hj) tells us what physical systems are at 
graph distance one. The set of such bound states as we vary j is the neighborhood of the system i. This is the set of vertices j a 
free particle in i can hop to. The projector has thus the meaning that the hopping interaction must be local in the sense just 
defined. 

Now we see that the term H ex is not satisfactory because exchange interactions are possible between any pair of vertices, no 
matter their distance. So quanta that are far apart can be converted in an edge between two points that were very far just before 
the conversion. Moreover, also the conjugate process is problematic, because it can easily lead to a graph made of disconnected 
parts. We implement locality by allowing exchange processes only between points that are connected by some other short path 
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of length L. Note that this refers to the locality of the state |\f» 9rap ' 1 ) a t time t relative to the locality of the state at time t — 1. 
Consider again the projector Pij on the edge (i, j) being present. Its L— th power is given by 

p ii= E Pik.Pk^.-.Pk^. (17) 

fei,...,fej_i 

For every state | € 'Hgraph, we have that Pij\^) ^ if and only if there is at least another path of length L between i and j. 
We can now modify the term H ex as follows: 

H ex = k ]T (S t j 1'^ ® (fctfttjfl + , } g, . (18) 

(hi) 

In the extended 5— level system, the exchange term is modified as \(s + l)(modS'))(s(modS')|(j ; j) <g> (bibj) R and similarly for 
the hermitian conjugate. 

This was the final step that brings us to the total Hamiltonian for the model which is 

H = H Hnk + H v + H ex + H hop . (19) 

In the following, we consider the theory for L = 2, which is the strictest notion of locality for the exchange interaction one can 
implement. 

B. Discussion of the model 

We can summarize the model in the following way. All we have is matter, namely the value of a function where the indexes 
i label different physical systems. We have chosen /; to be the number of quanta of the i-th harmonic oscillator. The bound state 
of a particle in i and a particle in j has the physical effect that other particles in i and j can interact. When there is such a bound 
state, we say there is an edge between i and j. Then other particles at i and j can interact, for instance, they can hop from i to 
j. The collection of these edges, or bound states, defines a graph which we interpret as the coding of the spatial adjacency of 
the particles (in a discrete and relational fashion). The physical state of the many body system is the quantum superposition of 
configurations of the particles and of the edges. The system evolves unitarily, and particles can hop around along the edges. But 
the distribution of the particles also influences the edges because some particles at vertices i, j can be destroyed (if i and j are 
nearby in the graph) to form another edge, and therefore making i and j nearer. The new edge configurations then influence the 
motion of the particles and so on. We have a theory of matter interacting with space. The intention of the model is to study to 
what extend such dynamics captures aspects of the Einstein equation and whether it (or a later extension of such a spin system) 
can be considered as a precursor of the gravitational force. From the condensed matter point of view, this is a Hubbard model for 
hopping bosons, where the underlying graph of the Hubbard model is itself a quantum dynamical variable that depends on the 
motion of the bosons. In the spirit of General Relativity, the edges (space) tell the bosons (matter) where to go, and the bosons, 
by creating edges, tell the space how to curve. 

We note that, in this theory, all that interacts has a local interaction by definition. We defined locality using the notion of 
neighborhood given by the set of systems interacting with a given system [? ]. We also note that, due to quantum superpositions, 
matter and space can be entangled. For this reason, the dynamics of the matter alone is the ruled by a quantum open system, the 
evolution for the matter degrees of freedom is described no more by a unitary evolution operator but by a completely positive 
map. We can show that the entanglement increases with the curvature. To fix the ideas, let us start with a flat geometry 
represented by the square lattice as the natural discretization of a two dimensional real flat manifold. In this model, a flat 
geometry with low density of matter can be described by a square (or cubic) lattice with a low density of bosons. This means 
that a particle is most of the time alone in a region that is a square lattice. The model will not then allow interaction between 
the particle and the edges, and all that happens is a free walk on the graph. On the other hand, when we increase the degree of 
the vertices by adding more edges, we make interaction, and hence entanglement, between edges and particles possible. This 
corresponds to increasing the curvature. In a regime of very weak coupling, i<fc(7,fi, entanglement will be possible only 
in presence of extremely strong curvature. From the point of view of the dynamics of the quantum system, this means that the 
evolution for the matter is very close to be unitary when curvature is low, while very strong curvature makes the evolution for 
the particles non-unitary and there will be decoherence and dissipation with respect to the spatial degrees of freedom. 

How does the graph evolve in time in such a model? The quantum evolution is complex, and since the model is not exactly 
solvable, numerical study is constrained to very small systems. In the next section we simulate the system with 4 vertices and 
hard core bosons. 

We can gain some insight from the analysis of the classical model, regarding H as the classical energy for classical variables. 
Since we delete edges randomly and build new edges as the result of a random walk of the particles, and there is nothing in this 
model that favors some geometry instead of others, we do not expect to obtain more than random graphs in the limit of extremely 
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long times. Indeed, we can argue as follows. With the exception of a very small number or graphs, all the other graphs belong 
to the set of graphs in which one can - under the evolution of our model - reach a ring. In practice, this means that there is a 
configuration in which one deletes all the edges without disconnecting the graph, and obtains many particles. This means that, 
starting from a state with zero particles, N vertices and Lq edges the number of eligible edges for deletion is Lq — aN with some 
constant a of order 1. So, as long as Lq > aN, the dynamical equilibrium between the number of edges and particles is realized 
when the rate of conversion of edges into particles equals the one of conversion of particles into edges. Let r be the number 
of edges destroyed (and pairs of particles created). If we assume that the particles move much faster than the edges, they will 
always be eligible for creating a new edge. The rates will be then the same when Lq—t = t which implies that at the dynamical 
equilibrium half of the initial edges are destroyed. Consider the case of the complete graph Kn, with zero initial particles. The 
initial number of edges is Lq = N(N — l)/2, and at long times, the dynamical equilibrium is reached when r = Lq/2. A 
similar reasoning can be applied to r instance, consider a square lattice of N vertices with periodic boundary conditions. In the 
state without particles, this is an eigenstate of the Hamiltonian, and therefore its evolution is completely frozen. Nevertheless, if 
we add a pair of particles, then all the other configurations of the graph can be reached, including the ring immersed in a gas of 
many particles. It turns out that the dynamical equilibrium is reached when Lq/2 = N of the initial edges are destroyed. The 
equilibrium state is obtained by deleting the edges randomly, and thus we expect to obtain a random graph. In order to obtain 
more interesting stable geometries, one has to put other terms in the Hamiltonian, that involve more edges together, meaning that 
curvature has a dynamical importance. The rigorous treatment of the asymptotic evolution of such graphs requires an analysis 
in terms of Markov chains, and it is developed in the next section. 



C. Trapped surfaces 

Another feature of this model is that allows for a very good trapping of particles and light, Let us consider a configuration 
of the system in which we have a region S of the graph with Ng vertices, that is highly connected, in an almost complete way, 
and then connected to the rest of the graph that can be a square or cubic lattice. We assume to be in the large limit for the ratio 
r between the number of edges in S and the number of edges between connecting S to the rest of the graph. Initially there are 
no particles. The region of the graph with "flat" geometry, that is, the square lattice, is frozen. In the highly connected lump, 
edges will start converting in particles. We have two time scales. The time scale to reach the dynamical equilibrium between 
the processes of destruction/construction of edges, and the time scale for particles to escape from S. We assume that r is large 
enough that the dynamical equilibrium is reached well before S starts losing particles. As a matter of fact, the region S will 
behave as the complete graph. Half of the edges will convert in particles, and we will have a dynamical equilibrium between 
a graph with O(Ng) edges in a gas of O(iVj) particles. Other particles coming from outside, will get trapped inside S too, 
therefore increasing the number of edges and particles within it. We want now to show that not even light can escape from this 
region. First of all, we want to have light propagating in our model. This can be done by adding to the Hamiltonian the term 
of Wen's U(l) theory for emerging light BUI . In the phase where the couplings of the U(l) theory are small with respect the 
others, we can have electromagnetic waves traveling on the graph. Now consider the state with very low density of particles, 
and with a graph that is represented by the wave function ^> = &a <E> *&b- Here B is a set of tib nodes such that ris *C riA- 
The wavefunction a is chosen to be the one representing a 3D cubic lattice. The wavefunction *S?b is instead chosen to be a 
very high dimensional hypercubic lattice, with dimension D ~ ng. Then one has to knit carefully the nodes in A with some of 

those in B. We choose to knit them with just of them, representing their "surface" E#. To summarize, the wave-function in 
the region outside the surface represents a discrete version of an euclidean three dimensional space. Inside the surface, we still 
have an euclidean space, but of very high dimension. On such a lattice, which has a well defined geometry, the emergent light is 
obeys the laws of geometric optics. On the other hand, a general wave-function $ for the edges would not posses any definite 
geometric meaning and the Wen's model cannot even be defined. 

The speed of light c on this graph can be estimated using the Lieb-Robinson bounds |25 1 and one can prove that c has different 
values in the different mediums A, B and it is proportional the geometric dimension of the medium so that we have cb/ca ~ «s ■ 
The argument about the asymptotic states shows that the number of edges in B will always be of the order of ris- In the geometry 
we have chosen, it makes sense to speak about the angle of a ray of light with respect to the normal to the surface separating A 
and B. The Snell's law of optics will imply that the critical angle for total internal reflection is 

8 C = sin -1 — ~ sin -1 ns (20) 
ca 

and therefore the probability of an emerging ray of light from B is of the order of n^. We see that light (and matter) are trapped 
within the hypersurface in the large ub limit. Eventually, some light can come out, and since the outer graph, even though less 
dense, contains more vertices and edges, eventually the region of space B has to evaporate. The whole process is completely 
unitary. Nevertheless, the emitted quanta of light and matter are entangled with the spatial degrees of freedom (the edges) inside 
Tib- But when B has evaporated, there are still edges there, with a density not much different from those outside Therefore 
the final state of the emitted quanta can be still entangled with the degrees of freedom inside £ b even if the "black hole" has 
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FIG. 1. Simulation of the system Hi/ 2 for N = 4. The initial state is pi = where \ipi) = |0000111111), that is, there are no 

particles and all edges are present. The parameters for this simulation are U = fi = l,t = k = .1. In the figure are plotted the quantities 
{SfSr} (red line), T{t) (blue line), S(t) (black line), P is (t), A(i)/3 ( green line) as a function of time. Revivals of the expectation value 
of the link operator coincide with revivals in the fidelity with the initial state. The initial value of the entanglement is 5(0) = because the 
initial state is separable. Notice that even though the fidelity is J-(t) > 0.85, the state has a non negligible entanglement. 



evaporated. The spectrum of the emitted radiation obtained by tracing out the spins in will therefore be mixed, even though 
the whole process is unitary. The black hole information paradox can be described in terms of entanglement. Pairs of particles 
inside/outside the event horizon are created and these pairs are entangled. The density matrix of the particles outside the horizon 
is therefore mixed, because there is a classical mixture of the particles coming from having to trace out all the degrees of freedom 
inside the horizon to which we have no access. The problem is that when the black hole disappears, there is nothing for the 
particles to be entangled with, and the mixture becomes a paradox. In our model, the particles are entangled with the spatial 
degrees of freedom. The disappearance of the black hole just means that the spatial degrees of freedom acquire a particular 
configuration, but the particles are still entangled with them, as in Eq.(fT0|i. 



III. THE MODEL WITH HARD CORE BOSONS 



A. Setting of the model 



In this section, we study the model Eq. ( 19 1 when the particles are hard core bosons. In this model, only at most one particle is 
allowed per site and the model can be mapped onto a spin system. We are particularly interested in the entanglement dynamics 
of the system. We have performed a numerical simulation of the time evolution of the model described by Eq. ( |19) . Since we are 
interested also in describing the quantum correlations in the reduced density matrix, we have resorted to exact diagonalization. 
In this way, we are able to compute the entanglement of the matter degrees of freedom with respect to the spatial ones. Of course, 
the simulation of a full quantum system is heavily constrained by the exponential growth of the Hilbert space. In this work, we 
have resorted to the simulation of hard-core bosons: at most one particle is allowed at any site. Hardcore bosons creation and 
annihilation operators must thus satisfy the constraints 

{b\) 2 = {U) 2 =0 (21) 
{Mj} = l (22) 
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FIG. 2. The initial state of this numerical integration is \ip2) = |0000111101), the parameters are in the "insulator" phase, U = fj, = 1, 
t = k = 0.1. The temporal scale is in units of h, on a range of 10 4 seconds. The diagonalization of the full system has been performed 
by means of Householder reduction, a) Time evolution of (Di(t)), (D^(t)). The damping of the oscillations is a sign of thermalization. b) 
Expectation values (Pi2(t)), {P2i(t)). The latter observable is fhermalizing. c) Von Neumann Entropy Si(t) for the sites i = 1, 2. We see 
that the entanglement dynamics is split in two different bands. The two vertices are only distinguished by the initial degree, d) Entanglement 
evolution S(t) and overlap with the initial state J-(t). The damping of J-(t) is a clear sign of thermalization. The entanglement S(t) between 
particles and edges shows the entangling power of the system, e) Expectation value of the particle operators at two different sites i = 1, 4. f) 
Concurrence C(t) as a function of time of the particles on the site i — 2 with the edge (2, A)(blue). Again we notice a damping of oscillations. 
Instead, the concurrence between the site 1 and the link 5 (red) is identically zero. 
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FIG. 3. The initial state of this numerical integration is = |0000111101), the parameters are in the "superfluid" phase, U = fi = 1, 
t = k = 0.1. The temporal scale is in units of h, on a range of 10 4 seconds. The diagonalization of the full system has been performed 
by means of Householder reduction, a) Time evolution of (Di(t)), (D^t)).. The oscillations have constant amplitude and the system does 
not present signs of thermalization. b) Expectation values of the link operators P14, P24- There is no sign of thermalization. The two values 
belong to two different bands depending on the initial value of the operator, c) Von Neumann Entropy Si(t) for the sites i = 1,2. We see 
that the entanglement dynamics is split in two different bands. The two vertices are only distinguished by the initial degree and the splitting is 
more marked than in the "insulator" case. Compare the result with the higher overlap of the operators Dij(t). d) Entanglement evolution S(t) 
and overlap with the initial state J-(t). Again the plots show no signs of thermalization. The behavior of J-(t) implies very long recurrence 
times, e) Expectation value of the particle operators at two different sites i = 1,4. f) Time evolution of Concurrence C(t). In blue is plotted 
the Concurrence between the vertex i — 2 and the edge (2, 4) for the "superfluid" case. Unlike the insulator case, the behavior of C(t) does 
not show any sign of thermalization. The concurrence between the site 1 and the link 5 is identically zero, as in the "insulator" case. 
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With these constraints, the bosonic operators map into the SU (2) generators 

S]o St (23) 
S- (24) 

b\h^ Q - S'^ (25) 

The local Hilbert space of a site i for a hard core boson is therefore that of a spin one half: 'H^ ch ~ C 2 . After the projection 
onto the hard-core bosons subspace, the model becomes a purely spin 1/2 model. For a system with n sites, the Hilbert space 
for the particles is thus the 2 Ar -dimensional Hilbert space Utosons — ^iLiH^ ~ <^ 2 ^> N _ The Hilbert space for the spatial 
degrees of freedom is still the 2 Ar ' Ar_1 ^ 2 -dimensional Hilbert space 

Hgraph = 1)72 He- The total Hilbert space is thus 

the 2 Ar ( Ar+1 '/ 2 -dimensional Hilbert space 

N N(N-l)/2 

Hspins = Hbosons ® Hgraph — (^) H\ (^) %e (26) 

i=l e=l 

As a basis for Hspins we use the computational basis. The basis is thus {|£i, £jv(iv-i)/2)ii! Jjv)}> where the first N(N ~ 
l)/2 indices label the edges of the graph, and the remaining N indices label the vertices. Of course if., = 0, 1 for every 
fc. 

After the projection onto the hard core bosons space H sp i ns , the model Hamiltonian becomes thus the spin one-half Hamilto- 
nian (for pi uniform): 

N /i \ 

ffi/2 =-u e s h) - a* E u " 5 " - f E P ^ ® + ^ r5 /) 

Let us examine the model in some limits. When the exchange term is vanishing, k — 0, the model has particle number 
conservation 

[H 1/2 ,N] = 0, iV = E 6 ^ ( 28 > 

i 

and therefore it has a U(l) symmetry, corresponding to the local transformation at every site given by 

M->n e< * 6 ' t6, M> <^[0,2tt) (29) 

while the total system with k 7^ does not have particle conservation because particles can be created or destroyed by means 
of the exchange term with the edges. Moreover, the k = system is self dual at /i = under the transformation 6, — > b\. For 
every separable state of the form \ip) = \ii, iN(N-i)/2) ® |V ) )6oso«s> me system is just the usual Hubbard model on the graph 
specified by the basis state iN{N-i)/2)- In the limit of —U positive and very large, all the edges degrees of freedom are 

frozen in the |1) state. The model becomes a Bose-Hubbard model for hard-core bosons on a complete graph. 

It is a typical feature of the richness of the Hubbard model that summing the potential and kinetic term gives a model with 
an incredibly rich physics. Depending on the interplay between potential and kinetic terms, it can accommodate metal-insulator 
transitions, ferromagnetism and antiferromagnetism, superconductivity and other important phenomena. The richness of the 
model comes from the interplay between wave and particle properties. The hopping term describes degrees of freedom that 
behave as 'waves', whereas the potential term describes particles |[T9l . As it is well known, the model is not solvable in two 
dimensions. The present model is even more complicated by the fact that the graph itself is a quantum variable. It is therefore 
extremely difficult to extract results from such a model. The hopping term in t favors derealization of the bosons in the ground 
state, while the chemical potential fi is responsible for a finite value of the bosonic density p in the ground state given by 



|E<^>- C3Q) 



The strength of determines how many bosons are present in the ground state. For // > 0, a large value of fi determines p = 1, 
meaning that the ground state has a boson at every site, whereas for p, < 0, a large value of p means there are no bosons in the 
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ground state p — 0. In any case, there is no possibility for hopping and this situation describe what is called a Mott insulator. 
On the other hand, for k = and t > p the hopping dominates and the system is in a superfiuid phase. The non vanishing 
expectation value in the ground state is that of the average hopping amplitude per link 

i,3 

We expect this situation to hold even for the weakly interacting system t ^> k ^ 0. As in the Hubbard model, there should be a 
quantum phase transition between the Mott insulator and the superfiuid phase for a critical value of p/t. An extensive numerical 
simulation of the ground state properties of the model is necessary to understand if, for k ^ 0, such transition belongs to the 
same universality class or a different one. It would also be interesting to understand whether there is a Lieb-Mattis theorem for 
such a system, namely that there are gapless excitations in the thermodynamic limit for the system of spins one-half. 

It should also be evident that depending on the interplay between potential and kinetic energy, the ground state of the system 
is entangled in the bipartition edges-particles. Starting instead from some separable initial state, the unitary evolution induced 
by H 1 / 2 will entangle states initially separable. 

B. Numerical analysis 

We have analyzed several aspect of the dynamics of the system in two different situations. The "insulator" case is the one in 
which the potential energies are dominant over the kinetic terms: U = p, = l;t = k = 0.1. The second situation is when the 
kinetic terms are much stronger, the so called "superfiuid" case: U = p = 0.1; t = k = 1. We have studied numerically the 
entanglement dynamics of the model, using as figures of merit the (i) Entanglement between particles and edges expressed by the 
von Neumann entropy S(t) of the density matrix reduced to the system of the particles, (ii) the Entanglement per site j expressed 
by the von Neumann entropy Sj (t) of the density matrix reduced to the system of just one site, and ( Hi) the Concurrence C (t) 
between a pair of edges, or particles or the particle-edge pair. This expresses the entanglement between these two degrees of 
freedom alone. 

We have simulated the system described by Hi/ 2 with N — 4 sites, which is 2 10 dimensional. We have labeled the 
sites i = 1, ..,4 starting from the lower left corner of a square and going clockwise. The basis states for the system are 
\J1J2J3J4', 614612623634624613) with Jj, = 0, 1. By direct diagonalization of the Hamiltonian, we compute the time evo- 
lution operator U(t) — e~ lHt . Starting from an initial state p(0), the evolved state is p(t) — U(t)pW (i). The entanglement 
S(t) as a function of time between particles and edges is obtained by tracing out the spatial degrees of freedom, we obtain the 
reduced density matrix for the hard core bosons: phcb{t) — T^graphp{t)- The evolution for the subsystem is not unitary but de- 
scribed by a completely positive map. The entanglement is computed by means of the von Neumann entropy for the bipartition 
W = Hbosons ® Hgraph, so we have 

S(t) = -Tr {p hcb (t) log p hcb (t)) (32) 

The single-site entanglement sj (t) is instead obtained by tracing out all the degrees of freedom but the site j and then computing 
the von Neumann entropy of such reduced density matrix. Finally, the last figure of merit to describe the entanglement dynamics 
of the model is the two-spins concurrence C(t) defined in ||29l . We define the r(i) reduced system of any two spins in the model, 
i.e., an edge-edge pair, or an edge-particle pair or a particle-particle pair. The entanglement as function of time between the two 
members of the pair is given by 

C{r(t)) = max(0, a/Vl- - \/% - y/Xi), (33) 

where Aj's are the eigenvalues (in decreasing order Ai > A2 > A3 > A4 ) of the operator T(t)(a y £g> cr y )r* (t)(a y ® <r y ). 

There are other important quantity to understand the time evolution of the model. We have computed the expectation value of 
the particle number operator S^S~ = |1)(1| j at the site i, the link operator P^, and the vertex degree operator Di = Pik> 
whose expectation value gives the expected value for the number of edges connected to the vertex i. The last important quantity 
is the fidelity J-(t) := \(tp(0)\ip(t))\ of the state \ip(t)) with the initial state \ip(0)). This quantity gives a measure of how much 
the state at the time t is similar to the initial state. 

The simulations have been carried out using two initial states \ r 4 > 2)- The state is the basis state describing the 
complete graph K4 without particles: = |0000111111). In Fig. [T|is shown the result of the simulation using as initial 
state, and for the model where the on-site potential energy is bigger than the kinetic energies, that is, in the "insulator" phase: 
U, p > t, k. Due to the very high symmetry of the Hamiltonian in the initial subspace, the system is basically integrable and 
we can indeed see a short recurrence time. Due to the initial symmetry of the state and the fact that no more than one particle 
is allowed at every site, the system is very constrained and it is integrable. The entangling power of the Hamiltonian is elevated 
and despite the fact that the overlap with the initial state is very high, the entanglement is non negligible. The expectation value 
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of every link is the same because of symmetry. For such an initial state, there is no qualitative difference other than different 
time scales between the "insulator" and "superfiuid" case. 

The time evolution starting from a just less symmetric state is far richer. The state \ip2) = |0000111101) is the basis state 
describing the square with just one diagonal, and again no initial particles. As anticipated, we have studied the model for two 
sets of parameters. The case (a) is the "insulator" case with parameters (i = U = l; t = k = 0.1. The case (b), or "superfiuid" 
case has parameters p, = U = 0.1; t = k = 1. 

"Insulator" case (a). — In the graph, Fig. [2ja) are plotted the time evolutions of (Dx(t)), (Z?2(i)) which have initial values 
of (£>i(0)) = 3, (1)2(0)) = 2. The oscillations of these operators are damped as well and the system is thermalizing towards a 
state which represents an homogeneous graph. It is very remarkable to see the phenomenon of eigenstate thermalization in such 
a small system. Recently, there has been a revival in the study of how quantum systems react to a sudden quench in the context 
of equilibration phenomena in isolated quantum systems, and our results are showing indeed that for such an isolated quantum 
system, the reduced system can thermalize due to the entanglement dynamics f38l [39l . In Fig. |2](b) are plotted the expectation 
value of the link operators 2 , P24 as a function of time In the initial state \ip 2 ) we have (V^I-P^IV^) = 1, (1^2^24 IV'2) = 0. The 
evolution of (Pi2(i)) is almost periodic but we see that on the other hand the oscillations of (P2<i(t)) are damping and the system 
is thermalizing. The behavior of the P13 operator is complementary to P24 and at long times \im t ^oo(P24(t) — Pi3(i)) = 0. 
In Fig.[2](c) we plotted the entanglement per site measured by the von Neumann entropy Sj(t) as a function of time. The sites 
considered are again j = 1 and j = 4. The two quantities split in two separated bands. Naively, one would expect that the vertices 
with higher degree are more entangled, but it is not so. Comparing with Fig|2](a) we see that the degrees (D\ (t)), {D±(t)) cross 
several time and have same time average. Surprisingly, the vertex with consistent higher entanglement is the one that started 
with higher degree at time zero, when the system was in a separable state. The system has a memory of the initial state that 
is revealed in the entanglement dynamics. This means that there are some global conserved quantities that are not detected by 
any local observable, but are instead encoded in the entanglement entropy, which is a function of the global wavefunction. The 
thermalization process is also shown in the behavior of the fidelity T(t) that presents damped oscillations, see Fig.[2](d). The 
behavior of the von Neumann entropy S(t) in Fig. |5](d) shows that the reduced system of the particles is indeed evolving as 
an open quantum system. Though some of the observables are thermalizing, the entanglement dynamics does not show any 
damping. In Fig. 0(e) is shown the time evolution of the expectation values of the number operators Ni(t) = S^~S~ (t)at the 
vertices i = 1,4. The vertices are distinguished by the initial state of the graph, namely by their degree. The average number 
of particles per site is Ni(t) ~ 0.05. This means that particles are basically involved in virtual creation/annihilation processes 
through which the graph acquires a dynamics. The next graphs show that the entanglement dynamics is all but trivial. The 
following Fig. |2] (f) shows the time evolution of the Concurrence C(t) between the vertex i = 2 and the edge (2, 4) which is 
in the state |0) at the initial time. We see deaths and revivals of entanglement, and the damping of the oscillation signals again 
some thermalization. 

"Superfiuid" case (b). — In this case, the kinetic terms k, t are dominant over the potential terms U, fi. The dynamics of this 
model for these parameters is completely different. We do not have any sign of thermalization. The degree expectation values 
(Di(t)), (£>2(i)) have a similar oscillating behavior with an even higher overlap, see Fig.[3](a). The link operators Pi4,P24, 
shown in Fig.|3](b) oscillate with no damping and are almost completely overlapped. To such distinct behavior with respect to 
the insulator case, we find a very strong similarity in the behavior of the entanglement per site Sj(t) (see Fig. [5](c) where, again, 
and in a more pronounced way, there is a splitting in two bands depending on the initial state of the system, and not on the degree 
(or other interesting observables) of the system during the time evolution. As in the insulator case, the vertex that started off 
with a higher degree is constantly more entangled than the one that started off with a lower degree, even if in the initial state they 
are both separable states and during the evolution all the relevant observables overlap strongly and have same time averages. 
This phenomenon again reveals how the entanglement contains global information on the state of the system that is not revealed 
in the usual local observables one looks at. The entanglement between edges and particles has a similar behavior than in the 
"insulator" case, but it is an order of magnitude greater, which is consistent with the fact that now the terms that couple edges 
and particles (and thus create entanglement) are larger. The superfluidity is revealed also in the behavior of the fidelity J-(t) that 
shows no sign of thermalization with constant amplitude of oscillations, see Fig.[3](d). The average number of particles per site 
is now Ni(t) ~ 0.52 and is homogeneous (Fig. pi (e)). Particles are delocalized over the quantum graph with a non vanishing 
expectation value. The Concurrence C(t) in Fig.l3](f) confirms that there is no thermalization in the system. 

To conclude this section, we have studied the model Eq.( 19 1 in the case of hard-core bosons. The usual Hubbard model with 
hard-core bosons on a fixed graph presents two quantum phases at zero temperature. An insulator phase, when the potential 
energy of the electrons is dominant, and a superfiuid phase, when the kinetic energy is dominant. In our model, the graph 
interacts with the electrons and the graph degrees of freedom are themselves quantum spins that can be in a superposition. 
We have studied numerically the entanglement dynamics of the system with four vertices starting from a separable state. The 
evolution with insulator parameters shows typical signs of thermalization in some of the relevant observables. Moreover, the 
entanglement dynamics reveals a memory of the initial state that is not captured in the observables. The behavior of the dynamics 
of the "superfiuid" system is completely different, in the fact that there is no apparent thermalization. The memory effect revealed 
by the entanglement dynamics is present in an even more pronounced way. There are many open questions to be answered: how 
the entanglement spectrum behaves and what it reveals of the system, what is the phase diagram of the model at zero temperature 



14 



in the thermodynamic limit, and a systematic study of the correlation functions in the model. We barely started studying the 
features of this model that presents formidable difficulties, but that promises to be very rich. 



IV. MARKOV CHAINS ANALYSIS OF THE MODEL 



In this section, we develop a general method to describe the evolution of graphs. We regard Eq.( 19 1 as the Hamiltonian for 
a classical model and consider a configuration of the system with a fixed number of edges and particles. The sum of these 
two quantities is a constant of the evolution. Moreover, it is safe to assume that almost all edges can be potentially converted 
in particles. The reason is simple: fixing the number of vertices, every connected graph (up to isomorphism) can be obtained 
by deleting and adding edges that are part of triangles. With this approximation, we expect that at long times a dynamical 
equilibrium is established between particles and edges. When considering the classical model, we can disregard superpositions 
and look at the dynamics as a discrete-time process with characteristics described as follows. For simplicity, we focus on the 
complete graph Kjy — (V(Kn),E(Kn)), with set of vertices V(Kn) and set of edges E(K^). Note that this constraint is not 
necessary, since we can start from any graph containing a triangle. The process simply needs at least one triangle in order to run. 
The process, starting from time t = 0, can be interpreted as a probabilistic dynamics gradually transforming the complete graph 
into its connected spanning subgraphs. These are subgraphs on the same set of vertices. Methods from the theory of Markov 
chains appear to be good candidates to study such a dynamics. We identify with a "graph of graphs" the phase space representing 
all possible states of the system considered. In this way a random walk on the graph, driven by appropriate probabilities, will 
allow us to study the behavior of the Hamiltonian, at least restricting ourselves to the classical case. Thus, the Hamiltonian 
transforms graphs into graphs. The Markov chain method suggests different levels of analysis: a level concerning the support 
of the dynamics; a level concerning the distribution of particles. At the first level, we are interested in studying graph theoretic 
properties of the graphs/objects obtained during the evolution, if we disregard the movement of the particles. This consists of 
studying expected properties of the graphs obtained by running the dynamics long enough. It is important to remark that the 
presence of particles does not modify the phase space, which, if we start from Kn, is the set of all connected graphs. The graphs 
obtained cannot have more edges than the initial one. Notably, particles only alter the probability of hopping between elements 
of the phase space. At the second level, we are interested in studying how the particles are going to be distributed on the vertices 
of the single graphs, and therefore in what measure the particles determine changes on the graph structure and consequently 
modify the support of the dynamics. When considering only the support, the Hamiltonian for the classical model determines the 
next process: 



• At time step t = 0, we delete a random edge of Go = Kjy and obtain G\. 

• At each time step t > 1, we perform one of the following two operations on G t : 



- Destroy a triangle: We delete an edge randomly distributed over all edges in triangles of Gt- A triangle is a triple of 
vertices {i,j, k} together with the edges {i, j}, {i, j}, {j, k}. 

- Create a triangle: We add an edge randomly distributed over all pairs {i,j} E(Gi) such that {i, k}, {j, k} £ 
E(Gt) for some vertex k. 

This process is equivalent to a random walk on a graph Qm whose vertices are all connected graphs. Each step is determined 
by the above conditions. The Hamiltonian gives a set of rules determining the hopping probability of the walk. The theory of 
random walks on graphs is a well established area of research with fundamental applications now ranging in virtually every area 
of science |28|. The main questions to ask when studying a random walk consist of determining the stationary distribution of 
the walk and estimating temporal parameters like the number of steps required for the walk to reach stationarity. The stationary 
distribution at a given vertex is intuitively related to the amount of time a random walker spends visiting that vertex. In our 
setting, the walker is a classical object in a phase space consisting of all connected graphs with the same number of vertices. 
Fig. |4]is a drawing of Q4, the configuration space of all connected graphs on four vertices. This is a graph whose vertices are 
also graphs. The initial position of the walker is the vertex corresponding to K4. 

The graph Qm (N > 2) is connected and bipartite. The number of vertices of C?jv equals the number of connected labelled 
graphs on N vertices. We need labels on the vertices to distinguish between isomorphic graphs. From the adjacency matrix of a 
graph G, we can construct the transition matrix T(G) inducing a simple random walk on G: [T(G)]ij — 1 /d(i) if {i, j} £ E(G) 
and [T(G)]ij = 0, otherwise. Here, d(i) :— \{j : {i,j} £ E(G)}\ is the degree of a vertex i. Notice that the degrees of the 
vertices in Qm are not uniform, or, in other words, Qn is not a regular graph. In fact, the degree of the vertex corresponding to 
Kn, which is the number of edges in this graph, is much higher than the degree of the graphs without triangles. In Fig. |4j it is 
easy to see that K4 has degree 6 and that the path on 4 vertices, drawn in the bottom-right corner of the figure, has only degree 
2. 

The evolution of a random walk is determined by applying the transition matrix to vectors labeled by the vertices encoding a 
probability distribution on the graph. The law (T(G) T )* Vq = v t gives a distribution on V(G) at time t, with the walk starting 
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FIG. 4. The graph <?4- The number of vertices is 38 and 72 edges. The number of vertices of Q n is exactly the number d n of connected labeled 

graphs on n vertices. The number d n satisfies the recurrence n2 |26|. The walker starts from the vertex corresponding 

to the graph K4. Even when we add particles, the graph Q n remains the support of the dynamics. The vertices of Q n are the possible states of 
classical evolution. In the quantum evolution, we have a weighted superposition of vertices. 

fromavertexi. The vector vj^ is an element of the standard basis of M. N . Thevectorv t = (Vj ,Vj , ...,Vj^) T is a probability 
distribution, being the probability that the walker hits vertex i at time t. The distribution 7r = (d(i)/2 \E(G)\ : i <E V(G)) 
is the stationary distribution, that is T(G)tt=tt. If G is connected and nonbipartite then lim^oo (T(G) T )* Vq = n l28l . The 
stationary distribution is independent of the initial vertex. Therefore, a walk on Qn can start from any vertex and the asymptotic 
dynamics remains the same. It is simple to see that Qn is bipartite. Then a random walk does not converge to a stationary 
distribution, but it oscillates between two distributions with support on the graphs with an odd and an even number of edges, 
respectively. In fact, for a bipartite graph G with V(G) = A U B, we have the following: lim t _ s . 00i even (T(G) T Y Vq 4 ' = Tr even 

with [Keyenli = d(i) / \E(G)\ if i € A and [ir even ]i = 0, otherwise; analogously for limt-^odd {T(G) T )\ { ^ = n dd- For 
instance, it follows that the stationary distribution of a random walk starting from any vertex of Q4 oscillates between the two 
distributions 




6 15 4 



and 

Tr even = (1/12, 0, -,0,5/72, .. .,5/72 ,1/31, .. .,1/31 ,5/72, .. .,5/72 ,0, 0, 0, 0). 

1 6 4 3 8 4 12 

In particular, for n odd , J2 i :d(i)=6l 7T oddli = \, ^t-.d^ii^oddh = g, ^-.d^^i^oddli = |; fonr even , J2i:d( i )=6i 7: odd\i = ^, 
E i :d( l )^Kddh = |, ■d(i)=2 [ L K odd\i — j2- Now, what is the most likely structure of a graph/vertex of Qn in which the 
random walker will spend a relatively large amount of time? In other words, where are we going to find the walker if we wait 
long enough and what are the typical characteristics of that graph or set of graphs? From the above description, one may answer 
this question by determining the stationary distribution of the walk in Qn- We do not have immediate access to this information, 
because we do not know the eigenstructure of T(Qn)- For this reason, we need some way to go around the problem. We can still 
obtain properties of the asymptotics by making use of standard tools of random walks analysis. In particular, as a first step, we 
are able to estimate the number of edges in the most likely graph in Qn- Even if this information is not particularly accurate and 
it is far from being sufficient to determine the graphs, it still can give an idea of their structure. The probability tt(G) that the 
walk will be at a given graph G after a large number of time steps is given, up to a small error term, by the stationary distribution 
7T of the walk. As we have mentioned above, n(G) is given by the number of possible transitions from G in the walk, divided by 
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a normalizing constant E which is independent of G. Based on this, it is possible to find the expected number of edges in a graph 
visited by the walk. We will provide here a sketch of the proof. A more extensive discussion is the Appendix. In total there 

are ) graphs with N vertices and k edges. By some of the classical results in the theory of random graphs ll27l we know 

that for k > ( 1+£ ^° s — ; the probability that a random graph is connected converges to 1 as TV grows. Let us recall briefly that a 
random graph is a graph whose edges are chosen with a fixed probability, equal and independent for each pair of vertices. The 
probability that a graph is visited by our walk will have k edges is given by J2v(g N )-\E(g N )\=k ^(Sn), where the sum is in fact 
taken over all connected graphs with k edges. In order to estimate this sum we must know how much it(Gn) can vary. However 
by the results mentioned earlier tt(Qn) cannot be larger than ( 2 ) / E and not smaller than (N — 1)/E, since that is the largest 
and smallest number of edges in a connected graph, and every transition in the walk can be associated with an edge in the current 

graph. Hence, the probability that the graph will have k edges will lie between (^) /E and ((?))(iV — l)/E. However 
for large N and k this value is completely dominated by the first term, which is of order 2 N /N for k = (^) /2. A more careful 
use of these estimates shows that the expected number of edges will be close to („) /2 ~ N 2 /A. Furthermore this results will 
hold true for any walk where the ratio between probabilities for the most and least likely graphs is not exponentially large in 
N, This observation tells us that the most likely graphs obtained during the process tend to have less edges than regular objects 
as lattice-like graphs. The number of edges in a square lattice with n 2 vertices is 2n (n + 1). For a cubic lattice on n 3 vertices 
this number is 3n (n + l) 2 , from the general formula dn (n + 1) ~ , where d is the dimension of the lattice. It seems natural to 
try to establish a relation between our walk and random graphs. After a first analysis, such a relation does not appear obvious. 
In fact, the walk on Qm is based upon a locality principle which is not usually defined when considering random graphs. An 
attempt to implement this principle for random graphs would consist in constructing Erdos-Renyi graphs starting from a random 
tree instead of the empty graph, that is, the graph with zero edges. A random tree, insures connectivity. Each pair of vertices 
at distance two is then joined with a probability p. If we keep adding and deleting edges, we obtain a dynamics similar to the 
one induced by our Hamiltonian. It is important to observe that the differences with the standard notion of random graph are 
essentially two: vertices at distance larger than two cannot be joined with a single step of the process; there is an additional 
probability of deleting edges. 

Let us keep in mind that so far we have not consider particles. Indeed, we have studied only a random walk on Q^, where 
this is the space of objects obtained by deleting and adding edges that form triangles. However, our Hamiltonian describes an 
evolution including particles. Each edge deletion creates two particles sitting at the end vertices of the deleted edge. These 
particles are free to move in the graph. Creation of another edge will depends on the number of particles. Only when two 
particles are located on two different vertices at distance two from each other, then we have a nonzero probability of creating an 
edge between such vertices and therefore creating a new triangle. Including particles, we can define the following process: 

• At time step t = 0, we delete a random edge of Go := Kn an d obtain G\. 

• At each time step t > 1, we perform one of the following two operations on G t : 

- Destroy a triangle: We delete an edge randomly distributed over all edges in triangles of Gt- When deleting an 
edge we create two (indistinguishable) particles. Each particle is located on a vertex of the graph according to the 
stationary distribution over Gt+i- This reflects the assumption that the particles thermalize. 

- Create a triangle: With a certain probability, we add an edge {i,j} ^ E(Gt) such that {i, k}, {j, k} £ E(Gt) for 
some vertex k. The probability of adding this edge is proportional to the probability of finding a particle at vertex 
i and a particle at vertex j at the same time t. When adding an edge, we destroy two particles. Specifically, the 
particles located in the two end vertices. 

Notice that the probability of deleting an edge is independent of the number of particles in the graph and their locations. On 
the other side, the probability of adding an edge is fundamentally connected to the number of particles. Higher is the number of 
particles in the graph G t E V(Qn) and higher is the probability of adding edges. The process exhibits a conservative behavior 
since the number of particles is always 

2^j- \E(G t )\. (34) 

So, the dynamics is again equivalent to a random walk on the graph Qn- This time the random walk is not a simple random 
walk, since the probability of each step is determined by the above conditions. In the transition matrix T(Gn) we can have 
[T(GN)]ij 7^ [T(GN)]i.k =h 0> whenever j ^ k. This fact gives different nontrivial weights on the edges of Qm- The transition 
matrix of the walk is then not necessarily symmetric and we need a normalization factor to keep it stochastic (i.e., the sums of 
the elements in each row is 1). Whenever an edge is deleted two particles are created. In the simplest version of our model all 
particles are distinguishable and at each time step of the walk all particles are redistributed according to a random walk on each 
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graph/vertex of Qm- By Eq. (34i, there are N 2p possible particle configurations. By the standard behavior of a random walk, 
particles will tend to cluster at vertices with high degree. For this model the state of the walk will consist both of the current 
graph G and the vector x of positions of all particles. Here the number of possible transitions will depend both on the structure 
of G, as before, and the number of particles. Since the number of particle configurations grows rapidly while the number of 
edges decreases, the walk will concentrate on connected graphs with few edges, rather than the denser graphs favored by the 
model without particles. If we make a rough estimate of the number of states corresponding to graphs with /2 edges, we 

see that they are fewer than 7y( 2)2(2) and that the number of states corresponding to graphs with O(N) edges are more than 

, for any e > 0. A comparison argument like the one used for the case without particles then shows that the expected 
number of edges for a graph/vertex will be o(N 2 ). It has to be remarked that an ad hoc tuning of the deletion probability for 
each edges should plausibly allow to obtain sparser or denser graphs. We have give a rough bound on the number of edges in 



a typical graph obtained via the Hamiltonian in Eq.( 19 1. The bound does not contradict the possibility that such a graph has 
an homogeneous structure like a lattice. Additionally to the number of edges, it may be worth to have some information about 
cliques. A clique is a complete subgraph. Cliques are then the densest regions in a graph. Knowing the size of the largest cliques 
gives a bound on the maximum degree and clearly tells about the possibility of having dense regions. In the Appendix we will 
prove that the growth of the largest cliques is logarithmic with respect to the number of vertices. This behavior also occurs for 
random graphs. 

Numerical simulations were performed to obtain information on the behavior of the classical system under different initial 
conditions. We are going to discuss the case of the complete graph Km as initial state. Complete graphs are interesting for 
several reasons. First of all, every edge of Km is eligible for interaction. This implies that edges rapidly transform in particles. 



As we can see in Fig. IV (a), the number of particles increases rapidly until it reaches an equilibrium value N${N). The number 
of steps to reach the equilibrium distribution is the same for all the graph sizes, and is of the order of the inverse of the only 
time scale introduced, given by ~ p~ . It is interesting to understand the equilibrium distribution of the degree for the various 
graphs Km- To obtain a better shape for this distribution, we increased the number of simulations from 30 to 60. The result can 
be seen in Fig. |lV](b). We find that the distribution is Poisson (D e is the degree), as it is for random graphs: 

1 (P e -f(iV)) 2 

P N {D e ) = -e 
it 



where R is a normalization constant. In Fig. IV (a) we find that the function f(N) is, for the graph Km, given by 

f(N) = f , 



while the function Q(N) of Fig. |rV|(b) is 

A' 



Q(N) . 

It has to be remarked that this result agrees for large N with the combinatorial proof in the last section. 



V. CONCLUSIONS 



The work presented here was motivated by the possibility that the answer to the problem of quantum gravity lies in the 
direction of emergent gravity, and, in particular, in a condensed matter perspective to emergent gravity. Along these lines, in 
previous work lfl6l [TTl . quantum graphity was proposed and analyzed, a background independent spin system for emergent 
locality, geometry and matter. The quantum states of this system are dynamical graphs whose connectivity represents locality. 
The Hamiltonian of [ 16 1 is not unitary: the universe starts at a high energy configuration (non-local) and evolves to a low energy 
one (local). This has clear limitations when applied to a cosmological context. The present model was originally intended as 
an energy-conserving version of 1161 . in which graph edges can be deleted, matter created and vice versa. The Hamiltonian we 
used is essentially an extension of the Hubbard model to a dynamical lattice. 

In our system, the basic building blocks of the theory are not events, but quantum physical systems Si, represented by a finite- 
dimensional Hilbert space Hi and a Hamiltonian Hi. That is, instead of an event, we have the space of all possible states of Si 
and the dynamical rules for the time evolution of these. The aim is to study the relationships among the Si and find geometry 
as the emerging structure imposed by these relationships, independently of their state. In fact, what we have can be considered 
as a unification of matter and geometry to matter only. Certain configurations (bound states) of matter play the role of spatial 
adjacency in the dynamical sense that they will determine whether interaction of other particles is allowed or not. In that sense, 
geometry is not fundamental but rather a convenient way to simplify the fundamental evolution of matter. From the condensed 
matter point of view, we have presented a Hubbard model in which the particles hop on a graph whose shape is determined by 
the motion of the particles and it is itself a quantum variable. 
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FIG. 5. a) Plot of the logarithm of the average number of particles N p as a function of the number of steps. The initial conditions were 
complete graphs with N vertices, Km, with hopping probability P h = 1 and interacting probability Pi = 0.1. b) In this plot we have the 
same quantity of (a), \og(N p (N)), plotted as a function of the graph size N at long times (when at equilibrium), c) Plot of the equilibrium 
degree distributions for each graph Kn. The hopping and interaction probabilities for each simulation are Ph = 1 and Pi — 0.1 respectively. 
The distributions are obtained averaging over 60 simulations. At equilibrium, we obtain Poisson distributions centered on y , as shown in Fig. 

|iv| (b), and variance Q(N) = ^f, as shown in Fig. Jivl (a). This a grees with standard results of the theory of random graphs. 
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This 



FIG. 6. a) Plot of the average values of the equilibrium degree distributions of Fig. i IV i (c). The red line is the fit, given f(N) = 
agrees perfectly with the various Kn- b) Plot of the variance of the equilibrium degree distributions of Fig. \TV\ (c). The red line is the fit, 
giving Q(N) = ^f. 



We simulated the quantum system for a complete graph with 4 vertices and hard core bosons and we analyzed the entanglement 
dynamics of the system, including the one between the particle and edge degrees of freedom. We argued that for the weakly 
interacting system, entanglement and loss of unitarity for the reduced system can be seen in presence of very high curvature. 
Moreover, the eigenstate thermalization of the model is studied under two different sets of parameters. Thermalization occurs 
when potential energy dominates over kinetic energy. 

The model discussed in this paper has some features of other models studied in a diversity of contexts and with different 
purposes. The first central aspect is a quantum dynamics involving a set of graphs. Quantum evolution on graphs is a subject 
of study related to spin systems, as a generalization of spin chains (see the review by Bose J33j), and quantum walks on graphs, 
where a particle undergoes a Schrodinger dynamics hopping between the vertices (see the review by Kempe [37]) where, at each 
time step, a particle is in a superposition of different vertices. Of course, the main difference with respect to our model is in 
the fact that in our model the graph is evolving in time as part of the system's evolution. Focusing on another area, it may be 
interesting to highlight a parallel with the work of Gudder [ 36 1 , who studied discrete-space time building on ideas of Bohm 1 32 1 . 
In the model described by Gudder, the graph is interpreted as a discrete phase-space in which the vertices represent discrete 
positions which a particle can occupy, and the edges represent discrete directions that a particle can propagate. This setting was 
primarily introduced to describe the internal dynamics of elementary particles. From this perspective, each particle is associated 
to a graph: vertices represent quarklike constituents of a particle and edges represent interaction paths for gluons which are 
emitted and absorbed by the vertices. Another aspect of our model is that we have many particles evolving at the same time. 
In the mathematical literature there is a growing number of examples of random walks with multiple particles/agents. A recent 
paper by Cooper et al. J34 1 studied properties of multiple random walks on a (fixed) graph assuming that the interaction between 
particles gives rise to various phenomena, like particles sticking together, annihilating each other, etc. It is important to mention 
that there are differences between walks with a single agent and walks with many agents. For instance, there are various settings 
in which k random walks visit all the nodes of a graph in expectation il(fc)-times faster than the case of a single walker [ 30 1 . 
Works addressing random walks on evolving graphs have been considered only recently, motivated by robotic exploration of the 
Web HQ [32. 

In a unitary system for cosmological evolution, the interesting question to ask is whether the system has long-lived metastable 
states. We studied this question rigorously, using analysis in terms of Markov chains, and also gave an intuitive argument, in both 
cases finding that dynamical equilibrium is reached when, starting from zero initial particles, half the initial number of edges 
are destroyed. In the model, there is a candidate configuration for a toy mechanism for attraction. A toy trapped surface can be 
modeled as the boundary of a region containing a highly connected subgraph. There is no singularity. To probe this region one 
can use Wen's hamiltonian for light EUll . with a coupling that is just a perturbation of the initial hamiltonian. We found that, in 
the limit of a large number of vertices inside the highly connected region, only a ray perfectly normal to the surface can escape. 

An important issue to emphasize before closing this paper is that our model assumes the existence of a notion of time and 
of time evolution as given by a Hamiltonian, as opposed to the constrained evolution of canonical pure gravity. It is a general 
question for all condensed matter approaches to quantum gravity whether such evolution is consistent with the diffeomorphism 
invariance of general relativity. While it is not possible to settle this question without first knowing whether the condensed matter 
microscopic system has a low energy phase which is general relativity, we can make a few comments, as well as point the reader 
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to more extensive discussion of this issue elsewhere ll22ll4D . In general, there are two possible notions of time: the time related 
to the goo component of the metric describing the geometry at low energy and the time parameter in the fundamental microscopic 
Hamiltonian. Let us call the first geometric time and the second fundamental time. In our geometrogenesis context, it is clear that 
the geometric time will only appear at low energy, when geometry appears. The problem of the emergence of geometric time 
is the same as the problem of the emergence of space, of geometry. The constrained evolution of general relativity, often called 
"time does not exist", refers to geometric time. By making the geometry not fundamental, we are able to make a distinction 
between the geometric and the fundamental time, which opens up the possibility that, while the geometric time is a symmetry, 
the fundamental time is real. It is important to note that the relation between geometric and fundamental time is non-trivial and 
that the existence of a fundamental time does not necessarily imply a preferred geometric time. We also note that, in the presence 
of matter in general relativity, a proper time can be identified. The particular system studied here has matter and in that sense it 
is perhaps more natural that it also has a straightforward notion of time. 

The toy model we presented here is very basic and there are several features we do not expect to see yet. For instance, there is 
nothing in the Hamiltonian to encourage the system to settle in metastable states that are regular graphs. Emergence of geometric 
symmetries such as Friedman-Robertson-Walker symmetries was not a goal at this stage but can be incorporated in future work 
by additional terms in the hamiltonian as in 11611 , or possibly also by introducing causality restrictions as in . Finally, we 
believe that this model is interesting from the condensed matter point of view. Condensed matter systems are always defined 
on a given lattice. In this model, the lattice itself is a quantum variable. Such a model, like the spin model of Eq. ( 27 1 can be 
realized experimentally in a system of quantum dots. We believe this is a novel way to think of condensed matter systems and, 
in perspective, potentially fruitful for the study of novel quantum phases. 

There are many potential generalizations and extensions of the model. Such extensions concern primarily the phase space and 
the nature of the particles. One possible generalization is to allow each edge to introduce a new vertex as a midpoint, possibly 
by absorbing particles. This gives homeomorphic graphs but with a growing phase space. The models could also be modified by 
allowing the degrees of the vertices to influence how the edges and the particles interact. A consequence would be a dynamical 
graph in which different regions are modified with different speeds. 
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VII. APPENDIX 



The following general result is useful to study the random walks described in the paper: 

Theorem VII.l Let G be a graph drawn from any random walk on Qn with stationary distribution it. Assume that there exists 
a < 2 such that 



Then 



1. E{\E{G)\) = \{ N 2 )+o{l) 



2. Fort> JN a 



q{tt) = max 

G,Hev(Q N ) [k] h 



< e N ° 



(35) 



Prob 



\E(G)\ 



1 fN 



2 V 2 



> t 



< 8exp(-Af Q ). 



Proof We will prove the theorem by proving a deviation bound for Prob 



\E(G)\-l( N 2 ) 



> t 



in terms of t. 



21 



Let 



At = <G 



\E(G)\- 



1 (N 



2 V 2 



< t 



Note that A t contains both connected and disconnected graphs. 

Let C be a function from the set of all graphs on N vertices to the set {0, 1}. Specifically, let 



C(G) := 



1 if G is connected; 
otherwise. 



Let C(N, t) denote the probability that a graph drawn uniformly at random from A t is connected. The probability that a graph 
drawn from the random walk belongs to A t is X^geA* G{G)[ir]c', likewise for G ^ A t . We will now bound the quotient between 
these two probabilities. It is well known, see e.g. [27 1, that 



for, e.g., t < N/2. Then 



C(N,t) = l + o(l) 



E G ^ t C(G)[7r] G - 9[7T) E G £A t C(G) 



(36) 



< 3(0 

<<?(0 



< 2ff(7T) 



\A t \C(N,t) 



This can be written as 



where 



23(0 



N,t 



N,t 



(37) 



Pjv.t = |A t |/2(") 

is the probability that a binomial random variable with distribution Bin((^) , |) deviates more than t from its expectation. Using 
the Chernoff bound we have then 



Pn,i < 2 exp 



-2t 2 
(a) 



If i 



7V Q (^) we can thus bound (|37|) as 



2sM ^— < 2o(7r)4exp(^^) 

= 8.g(7r)exp(-2n Q ) 
= 8exp(-iV Q ) 



(38) 



Since the denominator in the first step of Eq. ( |36[ l is less than 1 we have the following bound for our walk 

Prob [H £A t ]< G(H)[tt] h < 8exp(-7V a ) 



(39) 



H$A t 



For our range of a the value of t is o f which means that for a graph from A t the number of edges is | (^) + o(l) and the 
ected 

the total expectation is | (^) + o(l). 



contribution to the expected number of edges from graphs not in A t is between zero and (^)8exp(-7V Q ). This is o(l). Thus 
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Corollary VII.2 Let G be a graph drawn from a simple random walk on Qn. Then 

E(\E(G)\) = 1(^+0(1) 

Proof We know that the underlying graph of this random walk is bipartite and that in the asymptotic limit the stationary distri- 
bution oscillates between ir even and ir dd depending on whether we have taken an even or an odd number of steps. However if 
we start a new walk by not changing the graph in the first time step with probability 1/2, the new stationary distribution will be 

7T = -(n e ven + ^odd)- (40) 

Since we are looking at graphs with N vertices, d(i) is at most ( ^) and not less than N — 2. Hence, 

, . N 2 - N 



2n - 4 

and the corollary follows from Theorem |VII.l| 

Note that here we only used the fact the walk converges towards a stationary distribution on Qn and that g(n) is bounded by 
a polynomial in N for this walk. The same result holds for any form of walk on Qm for which g(n) is not exponential in N. 

Let us recall that a random event happens asymptotically almost surely, or a.a.s, if the probability for the event is 1 — o(l). 
A clique in a graph is a subgraph isomorphic to the complete graph. The clique number of a graph G, denoted by u>(G), is the 
number of vertices of the largest clique in G. 

Theorem VII.3 Let G be a graph drawn from a simple random walk on Qfj. Then there exist constants C\ < c 2 such that a.a.s 
the clique number oj(G) satisfies 

c log(iV) - Cl < u(G) < c log(JV) + c 2 . 

Proof It is well known, see, e.g., Chapter 11 of ||27| . that the expected clique number of a uniform random graph with edge 
probability \ is 

^log(iV) =c logiV 

and that the following concentration bounds hold 

Prob [u(G) - c logiV > r] < N~ r , (41) 

r — 2 

Prob [c log TV - uj(G) > r] < 7V~ L2 ^ J (42) 
We can now proceed in the same way as in the proof of Theorem |VII. l| using the set 

B? = {G\u(G)-cologN}\<t} 
to give a bound on the upper tail probability and 

Bi = {G\cologN-u(G)}\<t} 



for the lower tail probability, together with the bound on g(n) from the proof of Corollary VII.2 and the concentration bounds 



from the inequalities (41 1 and (42 1. This gives the following inequalities: 



N(N - 1) 

Prob KG) - c log N > r] < c x ^ _ - '- N~ r (43) 



Prob [c log N - u(G) > r] < ^ _ - > N L J (44) 

We can now use these inequalities to bound the contributions to the expected clique number. A clear but lengthly calculation 
can show that the contributions from the two tails are asymptotically bounded by two constants, giving us the bound stated in 
the theorem. 
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A rigorous analysis of the expected number of edges when we consider particles will be more difficult, since the model 
corresponds to a random walk on a directed graph, i.e., a graph in which edges have a direction. This is associated to an 
adjacency matrix which is not necessarily symmetric. There are transitions where, e.g., an edge is deleted and the particle 
distribution changes so that the endpoint of the edge do not have any particles on them, thus making the re-addition of the edge 
impossible in the next step. However, for states with a large number of particles, such as any graph with less than \ (") edges, 
the vast majority of transitions will be reversible, as it is in the case without particles. 

Let us consider the number of possible transitions from a state (G, x) defined as follows: here G has n vertices, t edges which 
are part of triangles, s edges with endpoints at distance 2, and p indistinguishable particles. If we assume that there are particles 
at all vertices we have three types of transitions which can be explicitly enumerated: 

1. There are transitions which correspond to redistributing the particles without changing G. 

2. There are transitions which correspond to deleting an edge and redistributing the p + 1 particles. 

3. There are s( n ~^^ 2 ) transitions which correspond to adding an edge and redistributing the p — 1 particles. 

Only the number of transitions of the last type is affected by the assumption that there are particles at all vertices. 

If p « c\ (™) we can estimate ( n+ £ _1 ) as 

{n/e) n c v " 

^ e l(l) n -\i + l + o{p- 2 )) 

n \nJ c 

Inserting this into the numbers of transitions given above, shows that for this model the maximum degree of the transition 
graph is bounded from above by a simple exponential, and the minimum degree is of course still greater than a multiple of n. 
According to Theorem 1 this is not sufficient to change the expected number of edges from being \ Q) . In order to make this 
analysis fully rigorous it is also necessary to show that the states with unoccupied vertices do not make a significant contribution, 
which will be lengthy but mostly a technical issue. 

The discussion for a case with distinguishable particles will be very similar but ( n+ p will be replaced by n p which is large 
enough to escape Theorem 1 . These counts also give an easy way to implement a simulation algorithm for both models. Just pick 
uniformly among all the possible transitions from the current state. If we do not consider particles or consider indistinguishable 
particles, the degrees of the vertices in the graphs will be close to those of random graphs with probability 1/2. This is more or 
less for the same reason that Theorem 1 works. The number of "typical" graphs in G(n, |) is so large that their behavior will 
still control these models. 
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